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Abstract. The roles of isospin asymmetry in nuclei and neutron stars are 
investigated using a range of potential and field-theoretical models of nucleonic 
matter. The parameters of these models are fixed by fitting the properties of 
homogeneous bulk matter and closed-shell nuclei. We discuss and unravel the 
causes of correlations among the neutron skin thickness in heavy nuclei, the 
pressure of beta-equilibrated matter at a density of 0.1 fm -3 , and the radii of 
moderate mass neutron stars. The influence of symmetry energy on observables 
in heavy-ion collisions is summarized. 
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1. The Symmetry Energy 

The nuclear symmetry energy is the energy required to create an asymmetry be- 
tween neutrons and protons (i.e., an isospin asymmetry). The preference for equal 
numbers of neutrons and protons is manifest in the difference between the deuteron, 
which is bound, and the dineutron, which is not. In nuclei, the preferred energy 
state would have equal numbers of neutrons and protons if the Coulomb interaction 
was absent. Expressed in terms of the energy density of homogeneous matter, the 
symmetry energy is given by 
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where n is the baryon density, 5 = (n n — n p )/n is the neutron-to-proton asymmetry, 
n n and n p are the neutron and proton densities, and e is the energy density. 

The symmetry energy influences several aspects of nuclear physics, from giant 
dipole resonances to heavy-ion collisions, and also several astrophysical processes, 
from supernovae to neutron stars. The broad influence of the symmetry energy is 
illustrated in Fig. Q][|I]. 
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T Fig. l r The multifaceted influence of the nuclear symmetry energy. , n 
In spire of trie diverse influences of tne symmetry energy, its> magnitude and 

density dependence are not well understood. Figure [21 shows the symmetry energy 
as a function of density for the models considered in this work (details are discussed 
in SectionyJ • The importance of the symmetry energy and our relative ignorance of 
it motivates studies of nuclear and astrophysical observables that relate to the sym- 
metry energy, and of how knowledge about the symmetry energy can help predict 
the outcome of experiments and astronomical observations. 



2. The Equation of State 

We will consider field-theoretical models in which a Walecka-type Lagrangian is 
utilized in the mean field approximation, and potential-models in which effective 
zero-range forces are used to construct an effective Hamiltonian density. In addi- 
tion to these two general classes, we also employ the equation of state of Akmal, 
Pandharipande, and Ravenhall [E] (APR). Because first principle Monte Carlo cal- 
culations of the structure of heavy nuclei are not yet available, we also perform 



Isospin Asymmetry 



3 




100 











0.1 



0.2 



0.3 



0.4 



0.5 



n (fm" 3 ) 



Fig. 2. The symmetry energy versus density for various equations of state. Solid 
(dashed) lines are for field-theoretical (potential) models. The thick solid line shows 
the APR symmetry energy for the low density phase. The arrow identifies the 
NRAPR (dotted line) and RAPR (solid line) models which have nearly identical 
symmetry energies at n = 0.5 fm -3 . 

field-theoretical (RAPR) and potential-model (NRAPR) fits to the APR EOS. De- 
tails of these fits are given in Ref. 

We require that all models reproduce the properties of nuclear matter at satu- 
ration density: 



In the case of potential models, calculations of nuclei are performed using the 
Hartree-Fock-Bogoliubov approach that includes pairing interactions. Hartree 
calculations [0] are employed in the field-theoretical approach, as a treatment of 
the exchange (Fock) terms is considerably more complicated than in the potential 
model approach. In both approaches, we require that the binding energy and the 
charge radii of closed-shell nuclei are reproduced to within 2% of the measured 
values. The scalar meson mass in the field theoretical approach was restricted to 
lie between 450 and 550 MeV. 



equilibrium binding energy : B 
equilibrium density : tiq 
incompressibility : K 
Landau effective mass : m* L 
symmetry energy : S v 



-16 ± 1 MeV, 
0.16 ±0.01 fm~ 3 , 
(200 - 300) MeV , 
(0.6 - 1.0) M , and 
(25 - 35) MeV . 
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In all cases considered the supranuclear EOS was constrained to yield a maxi- 
mum neutron star mass of at least 1.44Mq, the larger of the accurately measured 
neutron star masses in the double neutron star binary PSR1913+16 (see Ref. 
for a compilation of known masses). 

Monte Carlo calculations of low-density neutron matter have shown that the 
energy per baryon should be approximately half the Fermi energy [ E] . The APR 
equation of state, also based on Monte Carlo calculations, exhibits this behavior. 
By fitting the parameters of our field-theoretical model to APR, we can obtain this 
kind of behavior in RAPR as demonstrated in Fig. [3] where the energies per baryon 
of APR and RAPR are compared. 
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Fig. 3. The energy per baryon of neutron matter in the APR and RAPR mod- 
els. The energies in both models are approximately equal to half of the Fermi gas 
energies. 

3. Neutron Stars and Nuclei 

Recently, several empirical relationships have been discovered that underscore the 
role of isospin interactions in nuclei and neutron stars. 

The neutron skin thickness in nuclei is correlated to the pressure of pure neutron 
matter at sub-nuclear density : Typel and Brown [0 |Hj have noted that model calcu- 
lations of the difference between neutron and proton rms radii 5R — {r^) 1 ^ 2 — (r^) 
are linearly correlated with the pressure of pure neutron matter at a density char- 
acteristic of the mean density in the nuclear surface (e.g., 0.1 fm -3 ). We show this 
correlation in Fig.0]for the models considered in this work. The density dependence 
of the symmetry energy controls SR (termed the neutron skin thickness) in a heavy 
nucleus. 8R is proportional to a specific average of [1 — E sym (no)/ E sym (n)] in the 
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nuclear surface [ 1101 E| . To the extent that this correlation is valid, a mea- 
surement of SR will help establish an empirical calibration point for the pressure of 
neutron star matter at subnuclear densities. 
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Fig. 4. The neutron skin thickness SR of finite nuclei versus the pressure of f3- 
equilibrated matter at a density of 0.1 fm -3 . The blue squares display the results 
for potential models and the red circles show the results for field-theoretical modes. 
The circled triangles represent the potential (NRAPR) and field-theoretical (RAPR) 
model fits to APR. 

The neutron star radius R and the pressure P of neutron-star matter : Lattimcr 
and Prakash [El D2| found that the quantity i?P -1 / 4 is approximately constant, 
for a given neutron star mass, for a wide variety of equations of state when the 
pressure P of beta-equilibrated neutron-star matter is evaluated at a density in 
the range uq to 2no, where no denotes equilibrium nuclear matter density. Since 
the pressure of nearly pure neutron matter (a good approximation to neutron star 
matter) near ro is approximately given by n 2 dE sym /dn, the density dependence of 
the symmetry energy just above n determines the neutron star radius. Figure |3] 
shows this correlation as RP~ a versus R for stars of mass 1.4M for the EOS's 
considered here and densities n = 1.5 — 3no. In each case, the exponent a was 
determined by a least-squares analysis. 



4. Heavy-Ion Collisions 

4-1. Multi-fragmentation 

The breakup of excited nuclei into several smaller fragments during an intermed- 
iate-energy heavy-ion collision probes the phase diagram of nucleonic matter at 
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Fig. 5. The quantity RP~ a as a function of the radius R of a 1.4M Q star for 
pressures P determined at 3/2, 2 and 3 times equilibrium nuclear matter density. 
For each density, the best-fit value for the exponent a is as indicated. Circled stars 
indicate the results obtained with the APR equation of state. 

sub-saturation density and moderate (~ 10 — 20 MeV) temperatures. In this region 
of the phase diagram the system is mechanically unstable if (dP / dn)T. x < 0, and/or 
chemically unstable if (dfi p /dx)x,p < 0. (A pedagogical account of such instabilities 
can be found in Ref. [EI)- These instabilities, which are directly related to the 
symmetry energy at sub-saturation densities [114). are believed to trigger the onset 
of multifragmentation. Because of the instabilities, matter separates into coexist- 
ing liquid and gas phases, which each have different proton fractions, i.e. "isospin 
fractionation" [ 15 . This fractionation is observed in the isotopic yields which can 
potentially reveal information about the symmetry energy. Also, the scaling behav- 
ior of ratios of isotope yields measured in separate nuclear reactions, "isoscaling" , is 
sensitive to the symmetry energy [El 117) . To date, there are many suggestions of 
how the symmetry energy may affect multifragmentation [I18|. Ongoing research is 
concerned with an extraction of reliable constraints on the symmetry energy from 
the presently available experimental information. See also the work by Bao-An Li 
in this volume. 
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4-2. Collective Flow 

Nuclear collisions in the range Ei a f,/A = 0.5 — 2 GeV offer the possibility of pin- 
ning down the equation of state of matter above normal nuclear density (up to 
~ 2 to 3n ) from a study of matter, momentum, and energy flow of nucleons [I19|. 
The observables confronted with theoretical analyses include (i) the mean transverse 
momentum per nucleon (p x )/A versus rapidity y/y pro j [ I20| . (ii) flow angle from a 
sphericity analysis [[^, (hi) azimuthal distributions [E2> an d (iv) radial flow [I23|. 
Flow data to date are largely for protons (as detection of neutrons is more difficult) 
and for collisions of laboratory nuclei in which the isospin asymmetry is not large. 

The prospects of rare isotope accelerators (RIA's) that can collide highly neu- 
tron-rich nuclei has spurred further work to study a system of neutrons and protons 
at high neutron excess [EH ESI EE]< Observables that are expected to shed light on 
the influence of isospin asymmetry include neutron-proton differential flow and the 
ratio of free neutron to proton multiplicity as a function of transverse momentum at 
midrapidity. Experimental investigations of these signatures await the development 
of RIA's at GeV energies. In this connection, it will be important to detect neutrons 
in addition to protons. 

5. Outlook 

The correlations presented here will (with more progress in theory and experiemcnt) 
help to determine the symmetry energy and its density dependence. It will be 
interesting to explore further connections involving isospin asymmetry in heavy-ion 
observables and observables in nuclear structure and astrophysics. Work on this 
front is in progress [I27|. 
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